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Abstract 

We present the results of Monte Carlo simulations of the magnetic properties of a 
model for a single nanoparticle consisting in a ferromagnetic core surrounded by 
an antiferromagnetic shell. The simulations of hysteresis loops after cooling in a 
magnetic field display exchange bias effects. In order to understand the origin of 
the loop shifts, we have studied the thermal dependence of the shell and interface 
magnetizations under field cooling. These results, together with inspection of the 
snapshots of the configurations attained at low temperature, show the existence of a 
net magnetization at the interface which is responsible for the bias of the hysteresis 
loops. 
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1 Introduction 



In spite of its first report in the 50's [1], exchange bias effects in nanopar- 
ticle systems are still an active subject of research nowadays because of its 
technological application [2]. Current synthesis techniques have allowed the 
preparation of granular and particulate magnetic systems, in which ferromag- 
netic particles have either oxide shells or are coupled to AF environments. 
Shifts of the hysteresis loops after field cooling, similar to those observed in 
layered strucutures, have been experimentally observed in these systems [4]. 
However, despite the similarities in both cases, the theoretical models used to 
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account for this effect in layered systems [3] are not well suited for particle 
systems. Surface effects due to the reduced dimensionality of the nanoparticles 
and the exchange coupling at the FM-AFM interfaces are supposed to play 
a role in the observation of exchange bias that can be properly studied by 
Monte Carlo simulations of a single particle [5]. 



2 Model 

We have performed Monte Carlo simulations of a model of a nanoparticle 

— * 

consisting of Heisenberg spins Si on a simple cubic lattice. The particle shape 
is spherical with radius R = 12a (in units of the cell size a). The particle 
outer shell of spins having a width R Sh = 3a is considered as the surface and 
we will assign its spins magnetic properties different from those at the core. 
In particular, we are interested in a particle with a ferromagnetic (FM) core 
with uniaxial anisotropy kc along the long axis surrounded by a shell of pins 
with antiferromagnetic (AF) interactions and also uniaxial anisotropy ks- The 
interaction Hamiltonian considered can be written as: 

H /k B = -J2^s t -s J -k c J2(SD 2 

iec 

N 

-k s J2(s*y-Y,h-Si (f) 

iesh i=i 

where h is the magnetic field (h = fiH/ks will denote the field strength 
in temperature units, with \x the magnetic moment of the spin). The first 
term accounts for the nearest-neighbor exchange interactions, with different 
exchange constants. Core spin interactions are FM (J c > 0) whereas spins in 
the shell sorrounding the core have AF interactions (Jg < 0) with a reduced 
value with respect to that in the core in order to reproduce a Neel temperature 
lower than the Curie temperature of the core as in real core/shell oxidized 
nanoparticles [6]. Finally, in order to study the role of FM and AF interactions 
across the core/shell interface, for spins at the interface [i.e., the core (shell) 
spins having nearest neighbors in the shell (core)] the exchange constant is 

Jlnt^O. 



3 Results 

As a firt step towards the understanding of the exchange bias effects in mag- 
netic particles, we have started by simulating the field cooling process pre- 
vious to the hysteresis loop. For this purpose, we start from a disordered 
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configuration at a temperature higher than the Neel temperature of the an- 
tiferromagnetic shell and progressively reduce it in constant steps 5T = 0.1 
down to the measuring temperature of the hysteresis loop in the presence of 
a reduced magnetic field h — A K pointing along the positive z axis. At each 
temperature, the system is allowed to equilibrate during a number of Monte 
Carlo (MC) steps and after the magnetization is averaged over at least 10000 
MC steps using usual heat bath dynamics for continuous spins. The thermal 
dependence of the projection of the magnetization along the field direction is 
shown in Fig. 1, where the contributions of the spins in the core, in the shell 
and at the interface to the total magnetization have been plotted separately. 

We note that on field cooling the core magnetization tends to 1 indicating a 
progressive alignement of all the core spins towards the field direction driven 
by the FM coupling between them. In a similar way, the spins at the shell 
attain AFM order although with a small net magnetization originating from 
the uncompensated spin at the shell surfaces as can be more clearly seen in the 
snapshot presented in Fig. 2a. More interesting is the thermal dependence of 
the interface magnetization, which reaches a value of 0.37. This is an indication 
of the existence of uncompensated spins at the interface between the FM core 
and the AF shell of the particle. As it can be seen in the b-d panels of Fig. 2, 
the geometric structure and magnetic ordering of the interface in a core/shell 
nanoparticle is more intricate than in the case of FM/AF coupled bilayers with 
a roughness inherent to the geometry of the interface. For a nanoparticle, the 
core (shell) spins at the interface have different number of neighbors (1, 2 or 3 
shown in panels b, c, d of Fig. 2 respectively) on the shell (core) depending on 
their position. Moreover, the interface spins at the shell present regions with 
either local compensated or uncompensated magnetic order. 

Considering the configurations obtained after the previously described field 
cooling process as a starting point, we have also simulated hysteresis loops 
by cycling the magnetic field between the ±/i FC values at T = 0.1. Along the 
loops, the field is varied in constant steps 5h = 0.2 during which averages 
of the magnetization are performed during 100 MCS after other initial 100 
MCS are discarded for thermalization. The parameters that determine the 
characteristics of the hysteresis loops are the anisotropy constants at the core 
(which as been set to Kq = 1) and at the shell K^h, and the exchange constants 
at the core (Jc = 10, which fixes the Curie temperature Tc of the FM), at 
the shell Jsh and at interface Ji Ilt . It turns out that disorder and frustration 
at the surface induced by radial anisotropy and finite-size effects alone are not 
enough to produce sizable loops shifts as simulations performed for particles 
with no AF shell demonstrate [9]. Therefore, as most experiments in which 
exchange bias has been observed, the particle shell is in an oxidized state with 
Neel temperature below Tc, so that we have set Jsh = — 0.5Jc- Moreover, 
enhancement of the anisotropy at the surface as compared to that at the core 
is necessary to pin the shell spins into the cooling field direction so that a 
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shift of the loops occurs after field cooling, so that we have set K^h = 10), 
which also in agreement with the reported enhanced surface anisotropies due 
reduced local coordination at the outer particle shells [7,8]. 

As an illustrative case, we show in Fig. 3 the hystersis loops obtained by 
starting from a demagnetized state obtained after cooling the particle in zero 
field (ZFC) and after cooling in a magnetic field h^c — 4 (FC) for FM and 
AF coupling of the interface spins. As it is evident from the figures, the FC 
loops are shifted towards negative field values with respect to the ZFC loops 
and show a slightly increased coercivity independently of the sign of J\ nt . 
Changing the sign of the interface coupling affects the kind of magnetic order 
at the interface due to the direct coupling of the core to the shell through 
the exchange interaction (Jint), as is reflected by the different sign of the net 
interface magnetization just after the FC process (see lower panels in Fig. 3). 
For FM interface coupling (J Int = +0.5Jsh), the net interface magnetization 
after FC partially follows the reversal of the core. In contrast, for AF interface 
coupling (Ji n t = — 0.5Jsh), M\ nt remains always negative along the hysteresis 
loop as the magnetic field is not able to produce a positive net magnetization 
due to the strong AF coupling of the shell spins to the core. However, in both 
cases, the shift of the interface loops in the vertical direction indicates the 
presence of a fraction of spins that remain pinned during the reversal process 
and that are responsible for the observed shift of the hysteresis loops of the 
whole particle. 
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FIGURE CAPTIONS 



Fig. 1: (Color online) Thermal dependence of the magnetization of a core/shell 
particle when cooling form a disordered state at T > T N down to T = 0.1 in 
the presence of an external magnetic field h-pc = 4 K. The different curves 
correspond to the contributions of the core, shell and interface spins to the 
total magnetization Mpotai- 

Fig. 2: (Color online) (a) Spin configuration of an equatorial cut of the particle 
attained after the field cooling process described in Fig. 1. Core spins are dark 
blue, spins at the shell are green while inteface core and shell spins have been 
colored in light blue and yellow, (b-d) Configurations of the core (shell) spins 
at the interface having 1, 2 or 3 nearest neighbors in the shell (core). 

Fig. 3: (Color online) Main panels: hysteresis loops for a particle with R S h = 3a 
and radius R = 12a obtained from a demagnetized state (ZFC) and after 
cooling in a magnetic field hpc = 4. Lower panels: contribution of the interface 
spins of the shell to the total magnetization. The parameters used are: Jq = 10, 
J s = -0.5J C , K c = 1, K s = 10 and J Int = -0.5 (+0.5) in the left (right) 
panels. 
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